EQUILIBRIUM POINTS OF LOGARITHMIC POTENTIALS 
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Abstract. A notion of weighted multivariate majorization is defined as a 
preorder on sequences of vectors in Euclidean space induced by the Choquet 
ordering for atomic probability measures. We characterize this preorder both 
in terms of stochastic matrices and convex functions and use it to describe 
the distribution of equilibrium points of logarithmic potentials generated by 
discrete planar charge configurations. In the case of n positive charges we 
prove that the equilibrium points satisfy (j) weighted majorization relations 
and are uniquely determined by n — 1 such relations. It is further shown that 
the Hausdorff geometry of the equilibrium points and the charged particles is 
controlled by the weighted standard deviation of the latter. By using finite- 
rank perturbations of compact normal Hilbert space operators we establish 
similar relations for infinite charge distributions. We also discuss a hierarchy 
of weighted de Bruijn-Springer relations and inertia laws, the existence of 
zeros of Borel series with positive Z^-coefHcients, and an operator version of 
the Clunie-Eremenko-Rossi conjecture. 
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1. Introduction 

Let {zi}i^j\f be a sequence of distinct points in a domain f2 C C, where M Q N 
with \Af\ > 2. If at each Zi we place a particle with charge Ui - actually, an infinitely 
long thin rod perpendicular to the plane with uniform charge distribution - then 
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up to complex conjugation and a constant factor the resulting electrostatic force is 
given by 



(cf. \'2'M p. 10]). If \JV\ = oo then the series in 1)1. l|l and 1)1. 2|l converge absolutely 
for all z e f2 \ {zi}i^j^f provided that certain growth conditions are satisfied (see 
^3. HI below). In this case / is meromorphic while U is subharmonic in O. Since 
VJ7(z) = f{z) the critical points of the potential coincide with the equilibrium 
points of the field, that is, points upon which a free electron (or rod) once placed 
would remain. An important problem that goes back to Gauss is to describe the 
distribution and the relative geometry of the equilibrium points and charged parti- 
cles. One of the main purposes of this paper is to study this problem by means of 
the Hausdorff distance and the notion of weighted multivariate majorization that 
we define in ^|21as an extension of classical and matrix majorization. It is interest- 
ing to note that although our results are purely analytic or geometric in nature, 
they are obtained through a combination of the multilinear algebra and operator 
theoretic methods developed in [2111221 and |2SlESlEH|j respectively. The fact that 
such geometric properties could be established in this way was actually first noticed 
in loc. cit., where long-standing conjectures on the geometry of polynomials were 
solved by applying this type of techniques. Below we make extensive use of the key 
observation made in the aforementioned papers, namely that studying zeros and 
critical points of complex polynomials - or, more generally, poles and zeros of loga- 
rithmic potentials generated by finite positive charge configurations - is equivalent 
to studying spectra of normal matrices and their compressions. 

The weighted multivariate generalization of vector and matrix majorization in- 
troduced in a natural preorder on finite sequences of vectors in Euclidean space 
induced by the Choquet ordering for nonnegative Radon measures (Definition^. 
We characterize this new notion in terms of stochastic matrices and convex func- 
tions and we recover in this way several related results on classical and multivariate 
majorization. In the process we use a Hahn-Banach type argument to give a new 
simple proof of a fundamental theorem of Sherman (Theorem^. 

In f|31we state most of our main results. These describe the distribution and the 
relative location of equilibrium points of logarithmic potentials for discrete bounded 
planar configurations of positive charges and are of three different types. First, 
we consider n positive charges and prove that the equilibrium points satisfy (j) 
weighted multivariate majorization relations and are in fact uniquely determined 
by n — 1 such relations (Theorems j^lO • The higher de Bruijn-Springer relations 
thus obtained generalize some of the results of |25l 1^ to arbitrary (not necessarily 
equal) positive charges. Second, we show that the relative Hausdorff geometry of 
the extended set of equilibrium points and the set of charged particles is controlled 
by the weighted standard deviation of the latter ( Theorems QH^J. For equal charges 
the bounds given by these theorems are much better and also more natural than 
those predicted by the famous Sendov conjecture for complex polynomials and its 
converse, see, e.g., 0- As we explain in the sequel Q to this paper, these results 
also seem to suggest that the Hausdorff geometry is governed by natural disper- 
sion measures associated with the given charge configuration. Finally, we establish 
weighted multivariate majorization relations for essentially any bounded discrete 
set of positively charged particles with finite total charge for which equilibrium 




(1.1) 



and the logarithmic potential generated by these charged particles is 




(1.2) 
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points are known to exist (Theorems 0171). We further conjecture that equilibrium 
points actuaUy always exist under some mild assumptions (Conjectures EHSl in 

The proofs of the aforementioned results are given in 21 As we already pointed 
out, for finite charge configurations these are based on the multilinear algebra and 
matrix theory techniques of 25, 26, 28 . In the process we are led to consider several 
natural related questions, such as a hierarchy of inertia laws for finite systems of 
planar charged particles that we propose and discuss in fJJl (Conjecture This 
would yield a whole new series of weighted generalized de Bruijn- Springer relations. 
The techniques mentioned above actually extend to the case of infinite discrete 
distributions of charge, which allows us to make use of operator theoretic methods 
such as compressions of compact normal operators on separable Hilbert space in 
order to deal with this case. In this way we arrive at a great many related questions 
on the interplay between value distribution of meromorphic functions and variations 
of discrete spectra of compact normal Hilbert space operators under finite-rank 
perturbations. In ^HJEI give an operator theoretic approach to the well-known 
Clunie-Eremenko-Rossi conjecture JSI dealing with the existence of zeros of certain 
Borel series and we consider possible extensions of this conjecture as well as their 
potential implications (Conjectures EHEl- 

2. Weighted multivariate majorization 

The concept of majorization was first studied by economists early in the twentieth 
century as a means for altering the unevenness of distribution of wealth or income. 
The majorization preorder on n-tuples of real numbers - also known as the strong 
spectral order, vector majorization or classical majorization - was used by Hardy, 
Littlewood and Polya in their study of analytic inequalities ^Hl and by Schur in his 
work on the spectra of positive semidefinite Hermitian operators (cf. Classical 
majorization has since become an important tool in analysis, combinatorics and 
statistics and has recently found remarkable applications to quantum mechanics 
and entanglement transformations in quantum computation and information theory 
(see [2] and references therein). The strong spectral order essentially quantifies the 
intuitive notion that the components of a real n- vector x are "less spread out" than 
the components of another such vector y. Several matrix versions of this notion 
have been proposed and studied in various contexts (cf., e.g., [Hj and \'27\). 

In this section we generalize classical and matrix majorization by introducing 
a notion of weighted multivariate majorization that we shall use throughout both 
this paper and its sequel 0|. 



(2.1) 
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K- (2.2) 

n=l 

Given X = (xi, . . . , x^) e Xf^ we denote by X'^ = (xi, . . . , x^)^ its "transpose" , 
i.e., the ordered m-tuple consisting of the vectors Xi,...,Xm written in column 
form. We identify with the m x k real matrix whose i-th row consists of the 
coordinates of the vector x.^ in the standard basis of M.^. Let finally MJ^„ be the 
set of all row stochastic m x n matrices. 



4 



J. BORCEA 



The notion of weighted multivariate majorization defined below (see Definition^ 
is motivated by the following theorem of Sherman j tilj : 

Theorem 1. Let {X,a) e Xj'^ and {Y,h) e Xjj, where X = (xi, . . . ,x„i) e X^, 
a = (ai, . . . , Urn) G An, Y = (yi, . . . , y„) G A"^' and b = {bi, . . . , 5„) £ Ai- ^/le 
following conditions are equivalent: 

(i) For an^/ (continuous) convex function $ ; R'^ — > R one /las 



where X'^ is an m x k matrix and Y'^ is an n x k matrix obtained by some 
(and then any) ordering of the vectors in X'^ and Y^ , respectively. 



The special case when m = n and all weights are equal in Theorem ^ was long 
assumed to be an open question and appears only as an implicit conjecture in |27[ 
p. 433], which is the definite reference on majorization theory. Although essentially 
correct, Sherman's paper used rather complicated arguments and contained a series 
of misprints that required several subsequent errata, which may explain why the 
status of his proof seems to have remained unclear for such a long period. Below 
we present a short and simple proof of Theorem ^ 

Before embarking on the actual proof let us note that condition (i) in Sherman's 
theorem may also be written in terms of the Choquet ordering for positive measures, 
a notion that we shall actually use in ^ (see Definition 9 in loc. cit.). Indeed, any 
pair (X, a) G Xjj^ uniquely determines a discrete probability measure /^^^ on R*' 
concentrating a mass a.; at x^, 1 < i < m. Given two positive measures fii and 
/i2 on R'^ one says that /i2 dominates fii in the Choquet ordering or that /i2 is 
a dilation of fii, denoted fii -< /i2, if fJ-ii^) < M2('&) for any continuous convex 
function $ on K'°. Theorem ^(i) may therefore be rewritten as /^(^ a) ^ M(y-b)- 
We refer to ^1 ^1 ISHI for a detailed account of the Choquet ordering and 
various comparisons of nonnegative Radon measures and probability distributions 
in locally convex separable topological spaces. In particular, in the latter setting 
Sherman's theorem may in fact be viewed as a special case of the Cartier-Fell-Meyer 
theorem [3 123 13 and its subsequent extension due to Fischer-Holbrook . The 
most general result of this type, namely a basic representation theorem for convex 
domination of measures, seems to be contained in |13| . 

Proof of Theorem^ Since the implication (ii) ^ (i) is an immediate consequence 
of Jensen's inequality we shall only focus on the converse statement. Assume that 
condition (i) in Theorem 2] holds, let X = (xi, . . . ,Xm) G X^, a = (ai, . . . ,am) S 
•A-ni , Y (y 

ij • • • 7 Yri) G '^n, b = (fei, . . . , bn) G An, and define the following set 



Lemma 1. M{X,Y) is a non-empty closed convex subset o/MJ^„. 

Proof. It is clear that M {X, Y) is closed and convex. To show that M {X, Y) ^ % 
it is enough to prove that if Theorem ^ (i) holds then x^ e conv(yi, . . . ,y„) for 
1 < i < m, where conv(yi, . . . ,y„) is the convex hull of the vectors yi, . . . ,y„. 
Supposing that x; ^ conv(yi, . . . , y„) for some I G {!,..., to} let d denote the 
Euclidean distance in M'"' and consider the function $ : M''' ^ R given by ^(u) = 



m n 



^a.<f(x.) <^6,$(y^). 



(ii) There exists a matrix R G MJj^ „ such that 



X'^ ^ RY'^ and h = &R, 
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d(u, conv(yi, . . . ,y„)). Note that <i>(xi) > 0, ^{yj) = for 1 < j < n, and $ 
is a convex functfon since conv(yi, . . . ,y„) is a non-empty convex subset of R'^'. 
It foUows that J27=i'^i^i^i) X]j=i ^j*^(yj)' which contradicts condition (i) of 
Theorem^ Thus x^ e conv(yi, . . . ,y„) for 1 < i < m, as required. □ 

Given M = {rriij) e M{X,Y) we define a vector s(M) = (si(M), . . . , s„(Af)) e 
R" by Sj (M) = J27Li a^m^J, 1 < j < n. Set 

= {s(M) I Af e M{X,Y)}. 

Note that as(M) + (1 - a)s{N) = s(aM + (1 - a)iV) for any M, N e M{X, Y) 
and a E [0, 1], so that C{X, Y) is a convex subset of M". Moreover, since A4{X, Y) 
is closed the same must be true for £{X,Y). 

Lemma 2. With the above notations one has b G C{X,Y). In other words, there 
exists a matrix Mq € A4(X, Y) such that s{AIq) ~ b. 

Proof. Assume that b ^ C{X,Y). Since C{X,Y) is a closed convex set it follows 
from the Hahn-Banach theorem that there exist c G R and (ri, . . . , r„) e M" such 
that 

n n 

rjSj <c<Y^ rjbj for (si, . . . , s„) G C(X, Y). (2.3) 
Given an arbitrary vector y G conv(yi, . . . , y„) we define the simplex 



W{y)={ (Ai,...,A„)G [0,1]' 



Note that W{y) is compact for any y G conv(yi, . . . ,y„) since it is complete and 
totally bounded. We may therefore define a function G : conv(yi, . . . , y„) ^ R by 
setting 



G(y) = max ^r^ A^ 



(Ai, . . . , A„) G W{y) ) for y G conv(yi, . . . ,y„). 



Let Aj(y) G [0, 1], 1 < j < n, be such that G{y) = ^j^j(y)- Since Aj(yj) = (5^ 

one has G{yj) > r.j for 1 < j < n and thus 

n n 

Y.b,G{y,)>Y,h,r,>c. (2.4) 

i=i i=i 

On the other hand the left inequality in (|2.3I) implies that for any matrix M = 
{rriij) G M{X, Y) one has 

n m n 

c > Y^jSjiM) = ^^Oir^my. (2.5) 

Let us now consider the extension G of the function G to R'' given by G{y) = 
G{y) if y G conv(yi, . . . ,y„) and G(y) = for y G R'' \ conv(yi, . . . ,y„). It 
is not difficult to see that G is a concave function on R*^ . By Lemma one has 
Xi G conv(yi, . . . , yn), 1 < * < "t-, so that the numbers Aj(xi) are well defined for 
1 < I < TO, 1 < j < 71. Note that (Aj(xi)) G M{X, Y) and that by and 
the following holds 

m m n n n 

aiG(x,) = <^irj>^ji^t) < Y ^0^0 - ^J^(yi)- 
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This contradicts however condition (i) of Theorem ^ and thus we are done. It 
follows that b e F), which completes the proof of the lemma. □ 

The fact that (ii) (i) in Theorem ^ is equivalent to Lemma [3 □ 

Remark 1. We are grateful to the referee for pointing out the similarity between 
the Hahn-Banach type argument used above and the one given in the original proof 
of the Cartier- Fell-Meyer theorem |7]. 

Definition 1. The pair {X, a) e is said to be weightily majorized by the pair 
(y, b) e Xfj, denoted {X,a) -< {Y,h), if the conditions of Theorem ^ are satisfied. 

Remark 2. If a) -< {Y,h) then the a-barycenter of X must coincide with the 
b-barycenter of Y, that is, X^I^i ^j^i = Sj=i ^jyj- 

It is clear from Definition ^ that the weighted majorization relation is a preorder 
on X*^. Indeed, by Theorem ^ this relation is both reflexive and transitive. More- 
over, if TO = n and — bj — —, 1 < i,j < n, then Theorem ^ and Birkhoff's 
theorem [271 Theorem A. 2] imply that the weighted multivariate majorization re- 
lation induces a partial ordering on the set of unordered n-tuples of vectors in M.^. 
We shall refer to this partial ordering as the ordinary multivariate majorization 
relation, for which Sherman's theorem takes a particularly simple form: 

Corollary 1. If X ~ (xi,...,x„) e ^^^d. Y — (yi,...,y„) e X!^ then the 
following conditions are equivalent: 

(i) The inequality X]"=i — X]j=i*^(yj) holds for any convex function 
$ : M'^ ^ R. 

(ii) There exists a doubly stochastic nxn matrix S such that X'^ — SY'^ , where 
X'^ and Y^ are n x k matrices obtained by some ordering of the vectors in 
X'^ and Y^ , respectively. 

Remark 3. For fc = 1 Corollary ^ amounts to a well-known description of classical 
(vector) majorization due to Schur and to Hardy-Littlewood-Polya 123 • 

Remark 4. In the sequel we shall often consider weighted pairs of the form (^, a) 
and (y, b), where X and Y are unordered tuples of complex numbers. In this case 
we define the weighted majorization relation {X,a) -< (y, b), when appropriate, by 
identifying C with in Theorem ^ and Definition 

Another interesting consequence of Sherman's theorem is as follows: 

Corollary 2. Assume that m < n < k and let (xi, . . . , x„i) e X.}^^, (yi, . . . , y„) €E 
X^, (ai, . . . , a„) e An and (6i, . . . , 6„) G An- Form+l <i<n let (cji, . . . , Ci„) S 
An be such that b'^ := bj - J27=m.+i ^i^i] 7^ 0, 1 < j < n, and set Xj = Z]j=i CyTi- 
//yi, . . . ,y„ are linearly independent then the following conditions are equivalent: 

(i) Ya=i at'i'{yii) < bj^iyj) for all convex functions $ : R*' -> R. 

(ii) J2T=i aj*(xj) < &j*(yj) for all convex functions $ : R*' R. 
// either of these conditions is fulfilled then b'^ > for 1 < j < n. 

Proof. The implication (ii) (i) being trivial let us show that (i) (ii). Set 
X = (xi,...,x„) G X^, X' = (xi,...,x„), Y = (yi,...,y„), a = (ai,...,a„) 
and b = (6i,...,6„). Then condition (i) is equivalent to {X,a) -< {Y,h). Note 
that J2]=i^j — Y^^i'^i =: a e (0,1) and let a' = (a^^ai, . . . , a^^a™) G Am 
and b' = (a^^fo'^, . . . , q;""^6^). Using Theorcm^and the fact that yi,...,y„ are 
linearly independent one can easily show that b' G A^ and {X\ a') -< (Y, b'), which 
is equivalent to condition (ii). □ 
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Remark 5. For an arbitrarily given convex function $ : M —>■ R the inequality 
in Corollary |21 (i) is d priori weaker than the one in (ii) . As a special case let us 
consider three complex numbers 21,22,23 with 22^3 ^ R and a,/?, 7 e (0,1) such 
that 7 ^ (1 - a)/3 and 1 - 7 7^ (1 - q)(1 - /?). Then Corollary |3 implies that if 

a$(2i) + (1 - a)$(/322 + (1 - /3)^3) < 7^(22) + (1 - 7)^(^3) 

for all convex functions $ : C ^ M then 7 > (1 — a)(3, 1 — 7 > (1 — a){l — /?), and 
the inequality 

a$(2i) + (1 - a)(/3$(22) + (1 - PMzs j) < 7^(22) + (1 - jMzs) 
holds for any convex function $ : C ^ M. 

Remark 6. Theorem ^ generalizes several known results on vector and matrix ma- 
jorization, many of which seem in fact to be regularly rediscovered in various con- 
texts, see e.g. [0] and references therein. Two new interesting extensions of classical 
and multivariate majorization have recently been introduced and studied in |2()j . 



3. Geometry of equilibrium points: statement of main results 

In this section we formulate our main results. These describe the distribution of 
equilibrium points of logarithmic potentials for discrete configurations of positive 
charges (Theorems EllSl and Theorems EJIZIl and the Hausdorff geometry of the 
extended set of equilibrium points for finite charge configurations (Theorems El-El . 

3.1. Weighted majorization relations for finite charge configurations. Let 

Zi, I < i < n, he distinct points in the complex plane carrying positive charges a.i, 
1 < i < n, respectively, where n > 2. Throughout this section we assume wlog that 
the tti are normalized so that = 1- The logarithmic potential generated by 

these charged particles is given by 



[7(2) = ^ ai log 
1=1 

and the resulting electrostatic force is 

n 

f{z)^Y.^- (3-2) 



2 

1 

Zi 



(3.1) 



. Z- Zi 



The class of functions of the form H3.2|l is essentially the same as Sz-Nagy's class 
of generalized derivatives 32^. As explained in the critical points of U{z) coin- 
cide with the equilibrium points of the electrostatic field, that is, the zeros of f{z). 
By an argument reminiscent of the Gauss-Lucas theorem one can see that all equi- 
librium points lie in conv(2i, . . . , 2„), cf., e.g.. Proposition 3.1 in [26| . Remarkable 
generalizations of this simple fact were recently obtained in loc. cit., where it was 
shown that these points actually satisfy a whole sequence of majorization relations. 
Theorems I2H31 below further extend the aforementioned results. 

Notation 2. Let wi,. . . jWn-i be the critical points of U{z) counted with multi- 
plicity and recall the set Am defined in H2.1|l . For A € C \ {0} and ^ € C let 

W{\, ix) = {Xwi H- Ai, . . . , Xwn-i + m) e (3.3) 
Z(A, Ai) = (A21 -I- Ai, ■ • • , A2„ + Ai) e C", 

1 1 A , , / 1 — ai 1 — On \ 



1' 
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Theorem 2. There exists a matrix R S MJ^'Li „ such that W{\,^Y' — RZ{X,ii)^ 
and b = ai? for any X £ C \ {0} and G C, so that {W{X, ^),a) ^ (Z(A, /^), b). 
Equivalently, the inequality 

n—l n 

^$(Au;j+^) <^(l-a,)$(Az, +m) (3.4) 
is satisfied by any (continuous) convex function $ : C ^ M. 

Remark 7. The special case of H3.4|l when X — 1, ^ — and Oi — ^, 1 < i < n, is 
already a considerable improvement of the Gauss-Lucas theorem and was originally 
conjectured by de Bruijn and Springer ,6j. This conjecture was recently proved by 
S. Malamud |25l I2f)j and R. Pereira |2H1 independently of each other. 

Remark 8. Theorem|5]implies in particular that the following ordinary multivariate 
majorization relation holds in M.'^: 

(wi, . . . ,w„_i,C) -< (zi,...,0„_l,Z„), 

where C = J27=i '^i^i l*^!' ■ • ■ ' ari)-baryccnter of Z{1, 0) (cf. Remarks |21 and 0J. 

This relation was previously obtained in 26, Proposition 4.3] and 28, Theorem 5.4]. 
Note though that if n > 3 then even in the case of equal charges the above relation 
is in general weaker than the corresponding one given by Theorem |5] since the Zi 
are linearly dependent (compare with Corollary 

As we shall now see, (|3.4|l is actually but one relation in a sequence of (2) 
weighted multivariate majorization relations satisfied by the zeros and poles of 
f{z). Moreover, the equilibrium points of U{z) are in fact uniquely determined by 
n—l such relations. To formulate these results we need some new notation. 

Notation 3. Given fc G N and 1 < m < fc let Hk.m denote the mth elementary 
symmetric function on k symbols. For A G C \ {0}, G C, fc G {1, . . . , rt — 1} and 
1 < TO < fc set 

W^^{X,^i)^L^[r,,...,rk]{X,fi)) gC("^'), (3.5) 

where ^^[ri, . . . , rfc](A, ^) ^Ilk,m{Xwri + fi, . . . , Xwr^ +Ai), 

zW(A,A.) = (C™[si,...,sfe](A,M)) ec(^), 

where Cm[si, ■ ■ . ,Sfe](A, ^) = Ilk,miXzs^ + ■ • • , Az^^. + /i), 

»'"-((";T' 

Below we shall always assume that the components of Wm{X,^i), Zm (A, /i) and 
bl*^! are arranged increasingly with respect to the lexicographic order on the set of 
all indices consisting of fc-tuples of increasing numbers in {1, . . . , n}. In particular, 
for any 1 < to < (^) the TOth coordinates of Zm{X, /i) and bl*"'! are indexed by the 
same fc-tuple (si, . . . , Sk) satisfying 1 < si < . . . < < n. Note also that with this 
ordering (gSI) and (jHiSl imply that al^l = a, bl^l = b, M^f (A,^^) = W(X,fi) and 

z\'\x,^l)^zix,^i). 

Theorem 3.11 in [221 describes all possible sets of equilibrium points for loga- 
rithmic potentials generated by a fixed configuration of particles that are allowed 



G^(„-.), 



l<si <...<Sfc <n 
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to carry arbitrary positive charges. A natural related inverse problem is to describe 
geometrically the equilibrium points of a logarithmic potential associated to a given 
set of particles with prescribed positive charges. The following theorem solves this 
problem and provides a characterization of the zeros of the logarithmic potential 
U{z) defined in H3.1|l by means of n — 1 weighted majorization relations. 

Theorem 3. The following conditions are equivalent: 

(i) The zeros of f{z) counted with multiplicity are wi, . . . ,Wn-i, that is, 

n n 71—1 

^aiY\_{z ~ Zj) = - Wk), zeC. (3.6) 

1=1 j=i fc=i 

(ii) There exist A G C \ {0} and /i G C such that 

(M/fi(A,M),aW) ^ (zfi(A,Ai),bW), l<A:<n-l. 

(iii) There exist A G C \ {0} and fj, € <C such that for any k G {1, . . . , n — 1} one 
can find an ("^ ) x row stochastic matrix Rk{\, satisfying 

Wj!'\x,fxf = Rk{X,^i)zl^\x,^if and - aWi?fe(A,/i). 

(iv) There exist A G C \ {0} and /i G C such that for any k G {1, . . . , n — 1} the 
inequality 

J2 $(wfc[ri,...,rfc](A,Ai) 

l<ri<...<rfc<7i — 1 




< 

l<si<...<Sfc<n 

holds for all (continuous) convex functions $ : C — > R. 
(v) Either of conditions (ii)-(iv) is true for A = 1 and /i = 0. 

// the above conditions are fulfilled then for any k G — 1} there exists an 

("j^ ) X row stochastic matrix Rk such that 

M^f 1 (A, m)^ = (A, m)^ and bl^'l = aWi?^ (3.7) 

for all XeC\ {0} and ^ G C. 

Remark 9. Theorem extends Theo rem 4.7] to the case of arbitrary (not 
necessarily equal) charges, li at = ^, 1 < i < n, then Theorem |21 reduces to the 
characterization of the zeros and critical points of degree n complex polynomials 
by means of n — 1 majorization relations which was given in •26, Theorem 4.10]. 

From Theorem 13 we deduce a number of corollaries. 

Corollary 3. For any k G {1, . . . ,rt — 1} there exists an ("^^) x row stochastic 
matrix Rk such that 

W"^^ (A, ^if ^ RkZ^^ (A, iif and bl'^'l = a^ Rk (3.8) 

whenever 1 < m < k, X € C\ {0} and /i G C. In particular, for all such parameters 
k, m, X, /i and any convex function $ : C ^ M the following inequality holds: 



^ <^(um[ri,...,rk]iX,fi) 

l<ri<...<r^;<n — 1 

k 




(3.9) 

$(Cm[si,...,Sfc](A,M)). 
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Remark 10. Note that 1)3. 8() is d priori stronger than (|3.9() smce it asserts that 
there exist matrices Rk of the desired type that simultaneously realize the weighted 
majorization relations between the pairs ^W^m' (A, /i), a^*^!^ and ^Zm' (A, /i), b^*^!^ 
for all admissible parameters m, A and /i. 

Setting m = 1 in Corollary O yields the following generalization of Theorem |21 
which corresponds to the case when A; = 1 in H3.10|l below. 



Corollary 4. For any k G {1, . . . , n — 1} there exists an 
matrix Rk such that 

W[''\x,nf = RkZ[''\\,nf and b^ = aWi?fc 

for any A G C \ {0} and /i G C. Thus 



row stochastic 



iWj (zfkA,Ai),bWj , l<k<n~l, 
or equivalently, for all convex functions <f> : C ^ K and /cG {l,...,n — 1} one has 



(uJi[ri, . . . ,rfc](A,pi)^ 



l<si<...<Sfc<n 



(^Ci [si,...,Sfc](A,/^)y 



(3.10) 



Remark 11. The condition A may actually be omitted from the above inequa- 
lities. However, if A = the vectors defined in H3.5|l are independent of the Zi and 
Wj and the aforementioned inequalities become equalities. Thus the corresponding 
weighted majorization relations are trivially satisfied in this case. 

Taking A = 1, /i — and ^{z) — |2|" with a > 1 in Corollary we obtain the 
following moment inequalities for the charged particles and the equilibrium points. 

Corollary 5. For any k G {1, . . . ,n — I}, 1 < m < k, and a > 1 one has 

E \Ilk^,n{Wr,,...,Wr,)\°' 

l<ri<...<rfc<n — 1 

< |nfc,m(z5i,...,^sjr, 
l<si<...<Sfc<n 

where Tlk,m denotes as before the mth elementary symmetric function on k symbols. 
In particular, by setting Oi ~ —, 1 < i < n, and m = 1 in Corollary Owe get 



n — 1 
k 



E 



l<ri<...<r^. <n— 1 



< 



E 

l<si < . . .<Sfc <n 



(3.11) 



for all a > 1 and fc G {1, . . . , rt — 1}. Note that for k = 2 and a > 1 the inequalities 
in H3.11|l are quite similar to the inequalities 



n - 1 
2 



E 1^'- 

l<ri <r2 <n~l 



< 



that were conjectured by Rahman and Schmeisser in |3()l Remark 2.3.9]. In SjHlwe 
conjecture a whole hierarchy of weighted de Bruijn-Springer relations and inertia 
laws, see Conjecture ^ below. The latter actually contains as special cases all types 
of moment inequalities like those that we just described. 
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3.2. Relative Hausdorff geometry. The metric function between (closed) sub- 
sets of Euclidean space introduced by Hausdorff makes use of the following notions. 

Definition 2. The relative (alternatively directed or oriented) Hausdorff distance 
from a subset fli of C to another such subset Q2 is given by 

h{fli,il2) ~ sup dist(u', 5^2) = sup inf |w — zj. 

Other frequent names and notations are the forward Hausdorff distance from Qi to 
^2, hf(Qi,D,2) h(fli,n2), and the backward Hausdorff distance from Qi to fl2, 
hb{^i,i^2) ■— h{fl2,^i). The (symmetrized) Hausdorff distance between fli and 
O2 is then given by 

H{ni,n2) = max {h{ni,n2),h{n2,ni))- 

Let Z — {zi, . . . , Zn} be a set of charged particles and W — {wi, . . . , Wn-i} 
be the corresponding set of equilibrium points of the logarithmic potential defined 
in (|3.1|l . Theorems[21|2lyield a series of relations between the sets W and Z when all 
their elements are considered simultaneously. These relations go much beyond the 
obvious inclusion W C conv(Z) and suggest that no equilibrium point can actually 
lie too far away from the particles. Our next result gives a precise meaning to this 
intuitive fact and shows that the relative Hausdorff distance h{W, Z) is dominated 
by the a.-weighted standard deviation of the particles. The latter is given by the 
square root of the a.-weighted variance of z = (zi, . . . , z„), which we define as 

n n 

a2(z;a)2= inf V a,|z, - = V a,|2, - CP, (3.12) 

atEC ^ — ^ ^ — ^ 
i=l i=l 

where a = (ai, . . . , a„) g An and C stands as before for the a- weighted barycenter 
J27=i^i^i (cf. Remark IHJ. We shall actually prove a stronger statement involving 
the extended set of equilibrium points of the given system of positive charges, i.e., 
the set 

We^WU{C}^{wi,...,Wn-iX}- (3.13) 
Theorem 4. In the above notation one has h{W,Z) < h{yVe,2) < cr2(z;a). 

Remark 12. The bound given in Theorem 01 is sharp, as one can see for instance 
by placing n equal charges at the vertices of a regular n-gon. 

A natural question in this context is whether there is any analog of Theorem^ 
for the relative Hausdorff distance from Z to Wg. The following theorem provides 
an affirmative answer to this question when all particles are assumed to lie on the 
same line and shows that in this case the Hausdorff geometry of the set of charged 
particles and the extended set of equilibrium points is controlled by the weighted 
standard deviation of the particles. 

Theorem 5. // zi, . . . , z„ are collinear then H{Z, We) < cr2(z; a). 

Remark 13. Note that TheoremElis sharp and that a similar result cannot hold if 
We is replaced by W. To see this one may simply consider two distinct points in 
the plane carrying unequal positive charges. 

3.3. Infinite discrete distributions of charge. As we point out in fj^J several 
difhculties occur when trying to extend the results of i j3.1l to logarithmic potentials 
generated by an infinite discrete set of positive charges. To begin with, in this case 
already the fundamental question dealing with the existence of equilibrium points 
is quite delicate and, as a matter of fact, still an open problem (cf. |S1^^|^). 

Nevertheless, it turns out that in many situations one can obtain new informa- 
tion on the geometry of the electrostatic field through natural infinite-dimensional 
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versions of Theorem |5] As we shall now explain, this is for instance the case for 
collinear positive charges and, more generally, for essentially any bounded discrete 
set of positively charged particles with finite total charge for which equilibrium 
points are known to exist. Let {oijigN be a sequence of positive numbers and let 
{zijigN be a sequence of distinct complex numbers satisfying the conditions 

oo 

Oi > and \zi\ < p for i £ N, }^ at — 1, Hm Zi = p, (3.14) 

j=l 

where p is a fixed positive number. Using the same electrostatic interpretation as 
before - that is, placing a charge ai at Zi for i g N - we get the resulting force 



f[z)^Y.^ (3.15) 
and an associated logarithmic potential given by 

OQ 

U{z) ^^ai log 
1=1 

Set D(p) ^{zeC:\z\<p} and Df{p) ^ ©(p) \ {zJieN- By one has 



1 

Zi 



Eai 

so that the series in the right-hand side of H3.15|l converges absolutely for all z € 
D/(p) and thus / is a meroniorphic function in 0{p). One can see in similar fashion 
that U is subharmonic in 0{p). We have the following analog of Theorem [21 

Theorem 6. // {zijigN is a real sequence satisfying H3.14|l then the function f 
defined by (|3.15|l has an infinite discrete set of real zeros {uijjjgN md the sequences 
{zi}i£n and {wj}j£jq interlace on the real axis. Moreover, there exists an infinite 
matrix S = {sij)°°^i such that for all {i,j) G one has 

oo CXD OO 

Sij > 0, ^ Sik = 1, ^ sij = 1~ Oj, Wi = ^ SikZk- 

k=l 1=1 k=l 

In particular, for any A G C \ {0} and /i G C the inequality 

oo oo 

'^i^Wj +fi)< ^(1 - a,)<i>{Xzi + fi) 

j=i i=i 

holds for any convex function $ : C — > M. 

Our last main result is a natural analog of Theorem El for all bounded discrete 
configurations of positive charges for which equilibrium points are known to exist. 

Theorem 7. Let {zi}i^fi be a sequence of distinct complex numbers satisfying H3.14|l 
and let f be the meromorphic function in D(p) defined by H3.15|l . // / has a non- 
empty (discrete) set of zeros {wj}j,=fi then the inequality 

oo oo 

$(Awj +p)< ^(1 - a,)^{\z^ + n) 
]=i i=i 

holds for all X Cz C\ {0}, fi £ C and any nonnegative convex function $ : C — s- M. 

Remark 14. The fact that inequality (|4.18|) in M.3I mav be strict explains why the 
inequality in Theorem [71 holds in general only for nonnegative convex functions. 
The connection between majorization and such inequalities for nonnegative convex 
functions was first found in jl(i| . 
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Remark 15. Theoreni[71is reminiscent of Kato's bound for the variation of discrete 
spectra of Hilbert space selfadjoint operators under compact perturbations |21j . 

4. Proof of main results 

4.1. Proof of Theorems [gHU We shall use various methods from multilinear 
algebra and matrix analysis. These require a few notations and preliminary results 
and we refer to ^ and |27j for background material. 

Notation 4. Given m E N and 1 < fc < m let 

Qk,m = {i = {ii, ■ • ■ , ifc) I 1 < n < • ■ • < «fc < to}- 
If i? = (bij) is an TO X TO matrix and i,j G Qk.m we denote by 



Bii,j)=B 



«i, • ■ • 

VJlj • ■ • ijk/ 

the k X k submatrix of B lying in rows ii, . . . , ife and columns ji, . . . , jk- Below we 
shall always assume that the elements of Qk.m sue arranged in lexicographic order, 
i.e., if i and j are distinct fc-tuples in Qk,m then i > j if the first non-zero term in 
the sequence «i — ji, . . . , ifc — jk is positive. 

We fix an n-dimensional complex Hilbert space Ti. with unitarily invariant scalar 
product (•,•) and identity operator I = I-h = In ^ i(^), where L{'H) is the set 
of all linear operators on Ti.. Let A e LiTi) be a normal operator with spectrum 
Ti[A) = {zi, . . . , Zn} and choose an orthonormal basis (ei, . . . , e„) of Ti such that 
Aei = ZiBi for 1 < i < n. Define a unit vector 

n 

Vn = ^i^i (4.1) 

1=1 

and let P be the orthogonal projection on the subspace K. := v,^ of H. The operator 

A' = PAP\^ e L(/C) 

is called the P- compression of A jl|]). Recall the function / defined in (||-{.2|l . The 
following statement may be found in e.g. 156, formula (3.3)] and j2Hl Lemma 2.2]. 

Lemma 3. With the above notations one has 

{{A - z/„)- V„, v„) = ^''|^^r"/^j - -m forzeC\ S(A). 

Proof. Let (vi, . . . , v„_i) be an orthonormal basis of K.. The matrix representation 
of A' in the basis (vi, . . . , v„_i) is given by the (n — 1) x (n — 1) upper left-hand 
principal submatrix of the matrix representation of A in the orthonormal basis 
(vi, . . . , v„_i, v„) of Ti.. For any z g C \ S(vl) the {n,n) entry of the matrix 
representation of the (normal) operator [A — zln)^^ in the basis (vi, . . . , v„_i, v„) 
is given on the one hand by 

((A - z/«)-V„,v„) = ^^^e„f^(v„,e,)eA = -f{z). 

\,:=l ^ t=l I 

On the other hand, by Cramer's rule the (n, n) entry of the matrix representation 
of [A ~ zlfi)^^ in the basis (vi, . . . , v„_i, v„) of Ti is given by the cofactor of the 
(n, n) entry of the matrix representation of [A — zl-j-i) in the same basis. Thus 

for z e C \ S(y4), which proves the lemma. □ 
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Proof of Theorem\^ By Lemma 13 the zeros wi, . . . , Wn-i of / are the same as the 
eigenvalues of the P-compression A'. Let (vi, . . . , v„_i) be an orthonormal basis 
of K, that triangularizes A' . Then 

n 

^3 = v*A'v, = (A'v,,v,) - (Av,,v,) = ^z,|(vj,e,)|2 (4.2) 

1=1 

since = J^^^i (vj ! 6^)6^ for 1 < j < n — 1. Define the (n — 1) x n matrix R = (jij) 
by setting = |(vi,ej)p for 1 < i < n — 1 and 1 < j < n. Using the fact that 
(vi, . . . , v„_i, v„) is an orthonormal basis of TC one can easily check that 

n n 

E^»^- = E = ||v,||2 = 1, l<t<n-l, 

n— 1 n— 1 

E = E = ||e,||2 - |(e„ v„)P = 1 - a,, 1 < j < n, 

i=l 4=1 

SO that R e Mj;ii „ and b ai?. Note that by (|4.2() one has 

(Wi, . . . , W„_i)^ = . . . , Zn)"^. 

Hence (1^(1,0), a) -< (Z(l,0),b), which completes the proof of the theorem in the 
case A = 1 and /i = 0. The general case follows from this one by composing / with 
a non-singular affine transformation of the plane. □ 

Remark 16. The compression of a normal operator to a hyperplane is not necessarily 
normal. For instance, Fan and Pall jJS] showed that principal submatrices of a given 
irreducible normal matrix B are normal if and only if B has collinear eigenvalues. 

Remark 17. The proof of Theorem |21 given above is a slight modification of the 
corresponding arguments in 1211 ^^'^ EH]- 

We next prove Theorem |3| Clearly, conditions (ii), (iii) and (iv) in this theorem 
are equivalent by Definition and Theorem ^ Let us first show that (v) ^ (i). Set 

n — 1 n 

pi(z) = Y[(z - Zj), l<i<n, 
j=i i=i 

and note that condition (i) is equivalent to 

n 

g(»-i-'=)(0) = ^a,pi"-^-")(0), 1 < < n- 1. 
j=i 

Elementary computations show that 

k 

5(»-i-'=)(0) = (-l)'=(n-l-fc)! n^-.' 

l<ri<...<'r/s<n — 1 j — 1 
n / k \ k 

E (0) = {-l)\n - 1 - fc)! 1 - E n ' 

i—1 l<si <...<Sfc <n \ i—1 / i—1 

for 1 < A; < n — 1. Thus condition (i) amounts to saying that ifl<A;<n— 1 then 

k / k \ k 

E n-'-. = E i-E«- n-- (4-3) 

l<ri<...<rfc<n-l j = l l<si < . . .<Sfc <n \ i=l / i=l 

The fact that (v) (i) is now a consequence of the following lemma. 

Lemma 4. If condition (v) in Theorem\^ holds then (|4.3|l is true for I < k < n— 1. 
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Proof. Using e.g. condition (iii) in Theorem |21 with A = 1 and /i = one gets 

aWw^f 1(1,0)^ = b^zf 1(1,0)^, 1 < fc < n - 1, 
which is the same as the desired conclusion. □ 

We shah also need the following extension to normal matrices of a well-known 
result of Schur for Hermitian matrices (cf. QIEZI), see e.g. [23 Proposition 3.9 (a)]. 

Lemma 5. IfT — (tij) is a normal n x n matrix with eigenvalues ri, . . . ,t„ then 

{til, ■ ■ ■ , tnn) -< (ti, ■ ■ ■ , T„), 

i.e., {til, ■ ■ ■ , tnn)^ = *S'(''"i, . . . , Tyi)^ for some doubly stochastic n x n matrix S . 

Proof. Let f/ be a unitary matrix such that T — C/diag(Ti, . . . , Tn)U* and set S = 
U oU , where o denotes the Hadamard-Schur (entrywise) product of matrices and U 
is the complex conjugate of U . Then {tn, . . . , tnnY^ = S{ti, . . . , Tn)^ . It remains 
to note that 5' is a unitarily stochastic hence a doubly stochastic matrix. □ 

To complete the proof of Theorem 13 it will be convenient to work with certain 
suitably chosen matrix representations of the operators A and A' defined above. 
As in the proof of Theorem|211et us fix an orthonormal basis (vi, . . . , v„_i) of JC in 
which the matrix representation of the P-compression A' of A is upper triangular. 
For simplicity, denote the matrix representation of A in the basis (vi , . . . , v„_i , v„) 
again by A and let A^, 1 < i < n, be the degeneracy one principal submatrix of A 
obtained by deleting its i-th row and i-th column. Then clearly An = A' and so by 
LemmaOlthe spectrum of An is = {wi, . . . ,Wn-i\, where S(A„) is viewed 

as a multiset whose elements occur as many times as their algebraic multiplicities 
as eigenvalues of A„. Thus A is an n x n normal matrix with = {zi, . . . , 2„} 

such that 

c)' ^'-'^ 

where An is upper triangular with E(A„) = {wi, . . . , and C = {Avn, v„) = 

X]r=i '^i'^i ("^i' ■ ■ • ' an)-barycenter of Z{1, 0) (cf. Remark|HIl. Note that since 

the Wj appear on the diagonal of An Lemma IHl actually provides a new proof of 
Theorem|21 Let U — {uij) be the (unitary) change of basis matrix from (vi, . . . v„) 
to (ei, . . . , e„) so that we may write 

A ^ Udiag{zi, . . . , Zn)U* . (4.5) 

Proof of Theorem\^ Let us first show that if (|3.6|l holds and k E {1, . . . ,n — 1} 
then one can find an ("^ ) x (^') row stochastic matrix Rk such that 

M^f 1 (1, 0)^ = 1 (1, 0)^ and bl^'l = a^i?,.. (4.6) 

This is clear for fc = 1 by Theorem |21 and so we may assume that k >2. Let 

AW = (detA(iJ)),.,Q^^^ 

be the k-th compound matrix of A acting on the k-th Grassmann power a'^'H (recall 
that the elements of Qk n and therefore also the (^) x (^) determinants appearing in 
are arranged in lexicographic order). The Binet-Cauchy formula for compound 
matrices 27 Theorem 19. F. 2] and 1)4.5(1 imply that 

^C^) = {/Wdiag(zi, . . . , z„)W[/*W, 

which by an argument similar to the proof of Lemma [S] leads to 

{detA{i,i))lQ^^ = 5.zfl(l,0)^, Sk = C/C^-) o f/C^), (4.7) 
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where o denotes as before the Hadamard-Schur product. Note that Sk is an (^) x (^) 
doubly stochastic matrix since C/'^'^^ is unitary. Clearly, the diagonal elements of 
j^{k) corresponding to the k x k principal minors of A„ are of the form detA(i, i), 
i G Qfe,n-i- There are a total of ("^T^) such diagonal elements and these must 
actually coincide with the coordinates of the vector wj,*^' (1, 0)-^ since An is upper 
triangular. Let Rk be the ("■^^) x (^) row stochastic matrix obtained from Sk by 
deleting all elements S'fc(i, j) with i e Qk,n\Qk,n-i- We show that Rk satisfies H4.6|l . 
For this one has to prove that for any j = (ji, . . . ,jk) G Qk,n the identity 



E Sk 

l<'ii<...<«/^_i<n — 1 



holds, which is equivalent to 

E 

l<ii<...<ifc_i<n — 1 



detU 



or 



ii, . . . ,ik-i,n 



«i, ■ ■ ■ ,ik-i,n 
ji, • • ■ , jfc 

fc 



E' 



1=1 



E • 



1=1 



\detU{i,i)\' = l-J2 

i€Qk,n-l 



(4.8) 



(4.9) 



1=1 



By (|4.5|) one has 



/aj, 1 < J < so that by expanding each determinant 



with respect to the last row we get 

k 



detU 
and thus 
detU 



«i, ■ • • ,^fc-l,f^ 
ji,---,jk 



E(-i) 



'aj^ det U 



m=l 



«1, • • ■ ,»fc-l 
Jl J ■ ■ ■ J Jnn ■ ■ ■ 7 Jk 



«1, • • ■ , ik~i,n 
jl,---,jk 



E (-l)''^'\/«>airdetC/ 



r,s — 1 



«1, • ■ ■ , «fc-l 
Jl J ■ ■ ■ 7 Jr J ■ ■ ■ J Jfc 



detC/ 



Jlj ■ ■ ■ J Jsj ■ ■ ■ 1 Jk 



for 1 < ii < . . . < < n — 1. Therefore 



E 

l<ii<...<ifc_i<r! — 1 

where 



detU 



ll, . . . ,2A;-l,n 
jl,---,jk 



ar,s(j) = E 

l<ii<...<ifc_i<n— 1 



«1, • ■ • , «fc-l 
JIt ■ ■ ■ ■> Jrj ■ ■ ■ 1 Jk 



detU 



Jli ■ ■ ■ ■> J sj ■ ■ ■ 1 Jk 



Let us consider the fc x (n — 1) matrix 

V Jl7 ---Jk 

and the k x k matrix B = (bij) = AIM* as well as the k x k matrix given by 
its fc — 1 Grassmann power The Binet-Cauchy formula again 

implies that 

from which we deduce that 



^k+l-r,k+l-s 



ar,s{i), l<r,s<k. 
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Thus the left-hand side of (|4.8I) equals 

r,s— 1 

and by further manipulating this expression we obtain 

2 



E 

l<ii <...<ifc_i <n— 1 



detJ7 



Jl: 



■■■,3k 



^ detS(i,i) - fcdetS 



hence 

^ |detf7(ij)|^ = 1 - detB(i,i) + fcdetB. 

To show that the right-hand side of the above identity equals the right-hand side 



11 1 ■ 



a 



'3k I 



Since U is 



of 1)4.9(1 we proceed as follows. Define the vector u = {JH, 

unitary an easy computation shows that B = Ik — u^u, so that B depends only on 
the numbers aj^, I < I < k. It follows that the expression 

2 



E 



l<ii<...<j:fc<n-l 



detfj 



n, ■ 
Jl, ■ 



■ ,ik 

■ Jk 



= J2 Idet«7(ij) 



also depends exclusively on the numbers , 1 < / < fc, and is therefore independent 
of the Zi, 1 < i < n. On the other hand, by summing all coordinates corresponding 
to i e Qk,n-i in each side of ((4.7|l we get 



E 

l<ri<...<rfc<n-l j = l 




(4.10) 



E 



l<ii<...<ik<n—l 



detU 



,ik 

,Sk 



n 



Obviously, 1)3. 6|l defines the Wj, 1 < j < n — 1, as continuous functions of the 
Zi, ^ < i < n. By letting Zi for i G {1, ... ,71} \ {ji, . . . ,jk} and Zj^ 1 
for 1 < Z < fc in both (|CT))l and (g^J and comparing coefficients in the resulting 
identities we get (14.911 . which proves (|4.6(l . 

Finally, to prove the last statement of the theorem it is enough to consider the 
function f{X~^{z — fj,)), where A € C \ {0} and /i G C. The zeros and poles of this 
function are the eigenvalues of XAn + /i and XA + fi, respectively. These matrices 
satisfy relations similar to (|4.4() and (|4.5|l while the unitary matrix U is unchanged. 
One may therefore use the same arguments as above to get H3.7|l . □ 

Remark 18. Modulo some technical modifications that have to deal with unequal 
weights, the arguments given in the proof of Theorem|31are the same as those given 
by S. Malamud in his proof of jJEl Theorem 4.7]. 

Remark 19. An interesting byproduct of the proof of Theorem |3 is the following 
identity for a certain type of M-matrices. Let m S N, x = (.ti, . . . ,Xm) G (0, 1)™ 
with||x|p = J2^i^i < 1 and define the TO X TO matrix X = /m — x-^'x. Clearly, x'^x 
is a rank one matrix whose eigenvalues (with multiplicity m — 1) and ||x|p have 
modulus less than 1, so that X is an M-matrix (see, e.g.,|22]). By the arguments 
in the last part of the proof of Theorem Owe see that 



^ detX(i,i) - TOdetX = tr(/„i - X) 



both sides being equal to | |x| 
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Recall that if C is an n x n matrix and 1 < k < n one defines the k-th additive 
compound matrix Afe(C) as the (^) x (^') given by 



tC) 



(fc) 



t=o 



Equivalently, Afe(C) is the coefficient matrix of t in the expansion (/„ + tCY^'^ = 
/„ + tAfc(C) + 0(P). As an operator, this is the same as the restriction to a'^TY of 
the operator 

C «)/„«)••• (g) /„ + /„«) C (g) /„ (g) •••«)/„ + ... + /„ (g) /„«)••• (g) c 

defined on (g'^Ti. ft is well known that if C has eigenvalues Ai, . . . , A„ then the 
eigenvalues of Afc(C) are Xi^ + . . . + Xi^ , 1 < ii < . . . < ik < n (cf., e.g., P7]1. 

Proof of Corollary ^ Using (|4.5|l and arguing as in the first part of the proof of 
Theorem 13 we arrive at 



dt 



{In + tA) 



(fc) 



[/Wi. fdiag(l + tei,. 
dt \ 



A+tZr 



jj*(k) ^ 



Since An is upper triangular the diagonal elements of (/„ + tA)''^'^ corresponding to 
the k X k principal minors indexed by multi-indices i = (ii, . . . ,ik) G Qk,n~i have 
the form Jl^=i(l + twi^). The conclusion now follows from the obvious identities 



combined with the arguments in the second part of the proof of Theorem O 



□ 



Corollary ^ may alternatively be proved by using (|3.7|l with X = t and = 1 
and differentiating with respect to t. Using this alternative method we can actually 
prove Corollary 131 

Proof of Corollary \^ Let A G C \ {0} and G C. By (EHJl and Remark [TTl there 
exists an ("j^^) x (^) row stochastic matrix Rk such that the relations 

wl''\tX,tfi + lf = RkZl^\tX,tn + lf and = a^i^^ 
hold for all t e C. Thus 



^wi'^\tx,t, + if 



R,±-Zl^\tX,tf, + lf 



t=o 



for m e {f , . . . , fc}, which is the same as (jH.SII . 



□ 



Remark 20. The assumption that the points Zi, 1 < i < n, are distinct that was 
made at the beginning of H3.1l was not actually used in the above proofs. However, 
this condition is natural when viewing - as we did - these points as charged particles 
generating a resultant electrostatic force given by Clearly, the latter can be 

appropriately rewritten even in the case of multiple points. 

4.2. Proof of Theorems IIHSJ The inequality /i(W,Z) < h{yVe,Z) is trivially 
true since W C Wg. To prove the interesting part of Theorem 0] we make use of a 
result on spectral variations of normal matrices which is usually referred to as the 
Bauer-Fike theorem (cf., e.g., Theorem VI. 3. 3]). 

Lemma 6. Let A and B be n x n complex matrices. If A is normal then 

MS(S),S(A))<||A-B||, 
where h denotes the relative Hausdorff distance from S(-B) to S(y4.). 
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Recall the orthonormal basis (vi, . . . ,v„) of 7i where the matrix representation 
of the normal operator A with = Z is given by (|4.4(l . that is, 

In the above decomposition An is an (rt — 1) x (n— 1) upper triangular matrix with 
= W, C = (^v„, v„) = X;r=i and x,y e C"-\ Set B = A„ C- An 
easy computation now shows that 

n-l 

\\A~B\f = ||x|p = \\y\f = ^ |(Av„,v,)|2 = ||Av„||2 - |(Av„,v„)|2 

n 

= ^a,|z,|2-|Cp=a2(z;a)2, 

i=l 

which combined with Lemma yields the inequality stated in Theorem^ 

To prove Theorem \^ note first that by Theorem 21 it remains to check that 
h{Z,We) < C2(z;a). The latter inequality is invariant under non-singular affine 
transformations of the complex plane and it is therefore enough to show that 
h{Z,yVe) < cr2(z;a) whenever z = (zi,...,2„) e R". Since in this case A is 
a Hermitian operator we get from H4.11|l that x :— {xi, . . . ,Xn-iY — y and 
An = diag(wi, . . . , Wn-i)- Suppose Zi is such that 

min \zi — Wjl > (72(z;a) (4-12) 

l<j<n — 1 

and let u — (ui, . . . , u„)* S C" be a unit eigenvector of A with eigenvalue Zi. Then 
by (|4.11(l we may rewrite the identity Au = z^u as 

n-l 

— Zi)un — — and {wj — Zi)uj — XjUn, ^ < j < n — 1. (4-13) 

i=i 

Note that (72 (z; a) > since zi, . . . , z„ are assumed to be pairwise distinct. More- 
over, the fact that ||u|| = 1 combined with (|4.13(l implies that m„ ^ 0. From H4.12() 
and the last n — 1 identities in H4.13|l we get 

\uj\ < 0-2(2; a)"^|a;j||u„|, l<j<n-l, 

and then using the first relation in (|4.13|l we deduce that 

|C-z,| <a2(z;a)-i||x|p =a2(z;a). 

This shows that h{{zi}, We) < o'2(z; a), which completes the proof of Thcorcm[Sl 

4.3. Proof of Theorems |S|-[71 It is clear that up to a rescaling we may assume 
that p = 1 in both Theorem|Hland Theorem[71 Below we shall use the Hilbert space 
of square summable complex sequences Tl — P(N) with standard scalar product 
(•, •) and standard orthonormal basis {e^jigN. 

Proof of Theorem\^ We may assume wlog that {zijjgN is an increasing sequence. 
The fact that / has infinitely many zeros {lyj jjgN that interlace the poles {zi}i^j^ 
on the real axis follows by noticing that all zeros of / are necessarily real and 
that / is strictly decreasing on each interval (z^, z^+i), i G N. Alternatively, this 
may also be viewed as a consequence of the operator theoretic interpretation given 
below combined with Hochstadt's interlacing theorem for the eigenvalues of rank 
one perturbations of compact selfadjoint Hilbert space operators |2Uj . 

Let A denote the unique bounded operator on Ti. satisfying Aei = (1 — Zi)ei, 
I e N. Note that A is injective, selfadjoint and compact since the Zi are real and 
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Zi 1 as i OO. By assumption one has X^i^i ~ 1 we may define a unit 

vector V G 7i in the following way 



oo 



v = ^Va^'e^■ (4.14) 

i=l 

Set /C = V-'- and let P be the orthoprojection on JC. Recall from |22l §IV.6] that 
the Weinstein-Aronszajn determinant of the second kind associated with A and P 
is given by {{A — z/)~^v, v). A simple computation shows that 

((A-z/)-iv,v) zeC\{l-zJ,eN. (4.15) 

Let us now consider the compression A' oi A to JC, that is, the operator on /C given 
by A' = PAP\ic. Note that A' is itself a compact selfadjoint operator. We claim 
that the point spectrum of A' consists precisely of the zeros of /(I — z). 

Lemma 7. The set of eigenvalues of A' is {1 — WjjjgN. 

Proof. If u e /C is an eigenvector of A' with eigenvalue ( E C (actually R) then 

{A - C/)u = (/ - P)Au + {A' - C/)u {Au, v)v. (4.16) 

From (|4.16(l we see that {Au, v) 7^ since otherwise u would be an eigenvector 
of A, say u = for some i E N, which would then lead to the contradiction 
(Au, v) = {Aei,v) = (1 - Zi)^, ^ 0. By iflT^ and (f^^ we get 

/(I - C) = ((A - C/)-^v, v) = (Au,v)-i(u,v) = 

so that 1 — C = for some j E N. Conversely, let w be such that /(I — w) = 0, 
that is, ((A — u'/)~"'^v, v) = 0, and set u = (A — u;J)~^v. Then clearly u 7^ and 
one has 

V = (A - w/)u = {Au, v)v + (A' - wl)u. 
Since (A' — w/)u E JC and v e JC^ it follows that 

(A' - wl)u = (1 - (Au, v))v = 0. 
Thus u is an eigenvector of A' with eigenvalue w. □ 

It is not difficult to see that A' is injective. Indeed, if 7^ u g /C is such that 
A'u = then from PAu = and the fact that A is injective we get Au = av for 
some a € C \ {0}. Hence 

0= (u,v) =a(A-iv,v) =a/(l), 

which is a contradiction since the assumptions of the theorem (more specifically, 
condition (|3.14(l with p = 1) imply that /(I) > 0. 

By the spectral theorem for compact selfadjoint operators ^2 Corollary X.3.3.5] 
there exists a complete orthonormal system {f,}jgN in JC consisting of eigenvectors 
for A'. Set Sjk = |(f, ,efe)|2 for j,kE N. Then 

00 

1 - W, - (A'f„ f,-) = (Af„ f,) = ,5,fc(l - Zfe). 

fe=l 

Both {v} U {fjljgN and {efe}jgN are complete orthonormal systems in TL, so that 

00 00 

00 00 

Y^s.k = ^ |(f„efe)P = ||efe||2 - |(v,eOP = 1 - flfe, A; e N, 
which completes the proof of the theorem. □ 
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Proof of Theorem^ Let {zijigN be a sequence of distinct complex numbers satis- 
fying H3.14|l with p = 1 and let A be the injective compact normal operator on H 
defined by Ae.i = (1 — Zi)ei, i E N. As before we denote by A' the compression 
PAP\)c of A, where K. = v"*" and v is given by H4.14|l . Clearly, A' is a compact 
operator on JC. Moreover, since 3fi/(l) > the arguments in the proof of Theorem|Sl 
show that A' is injective. Note that Lemma |7| remains valid in the present setting 
and that by assumption / has a discrete set of zeros {wjjjgN- It follows that the 
point spectrum of A' is {1 — Wj}j^n- Denote by JCa' the closed linear hull of the 
root subspaces of A' and let A' be the operator on K-a' induced by A' . Recall the 
operator version of Schur's lemma given in 17, Lemma 4.1]. 

Lemma 8. With the above notations there exists an orthonormal basis {fjjjgN of 
K-A' in which the matrix of the operator A' has triangular form and (A'fj,fj) = 

1 - wj, j e N. 

Set Sjk — |(fj,efc)p for j, fc e N. By Lemma|Slwe have 

oo 

fc=l 

It is clear that 

oc oo 

^s,,=^|(f„efe)p = ||f,.|p = i, jen. (4.17) 

fc=l k=l 

On the other hand {v} U {fjjjgN is an orthonormal system in H so that by Bessel's 
inequality we get 

oc oo 

=^|(f,,efc)|2 < llefclp- |(v,efc)|2 = l-flfe, keN. (4.18) 
i=i i=i 

The inequality stated in the theorem by means of nonnegative convex functions is 
an immediate consequence of H4.17|l - H4.18|l . □ 



5. A HIERARCHY OF DE BRUIJN-SpRINGER RELATIONS AND INERTIA LAWS 

Theorems suggest even deeper physical connections between the set of equi- 
librium points and the given set of charged particles in the case of finite planar 
systems, such as a whole new hierarchy of inertia laws for planar solid bodies as- 
sociated with these sets. In order to make a precise statement let us identify the 
complex plane C with the xiX2-planc in M"^ equipped with the standard coordinate 
system (O; xi, X2, x^). For an integer fc > 2 we denote by the standard simplex 

Sfc = {{ti, . . . ,tfc) e [0, 1]^ I ii + . . . + tfc = 1} . 

Definition 3. Let i be a line in M^. Given a convex fc-gon /C in C with vertices 
vi, . . . ,Vk and a > 1 we define the a-moment of IC with respect to L to be 

Ia,LilC) ^ J'J d{tivi + ...+tkVk,L)°'dti---dtk, 

where d{-,L) denotes the distance function associated with L. 

Remark 21. If /C is viewed as a planar solid body with uniform mass distribution 
then Ia,L{K^) may be interpreted as a higher moment of inertia (or generalized 
angular momentum) of K. with respect to the rotation axis L. 
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Let Sm denote the symmetric group on m elements. It is not difHcult to see that 
if K, is as above and ui, . . . , Um S C are such that /C — conv(ui, . . . , Um) then 

d{tT,(l)Ui + . . . + L^dti ■ ■ ■ dtm = la.LQC) (5.1) 

for any tt g so that Ia,L{^) depends only on the extreme points of /C. 

The following conjecture is a powerful generalization of Theorems |2H3 and sug- 
gests a "hierarchy" of weighted de Bruijn-Springer relations and inertia laws for 
planar polygons constructed from the equilibrium points and the charged particles. 

Conjecture 1. Let Zi, 1 < i < n, be distinct points in C and > 0, 1 < * < n, be 
such that X]"=i — 1; where n>2. Denote by Wj, I < j < n — 1, the zeros of 



Z ^ y — 



z - 

i=l 



and let ^ : <C be an arbitrary convex function. Then for any k £ {l,...,n— 1}, 
1 < m < k, and (ti, . . . , tk) G C*^ one has 



l<ri<...<rfc<n-l TreS 




^ (n/£,m(t7r(l)^si, ■ • • , ^7r(/£)Zsfc)) : 



(5.2) 



where Ilk,m is the mth elementary symmetric function on k symbols. In particular, 
for any line L in R"^ and a > 1 one has 



Ia,L (conv(wri , . . . , J) 




(5.3) 

1 - y^ fls, laX (conv(zsi, . . . , ZsJ) • 




l<ri<...<rfc<n — 1 



The fact that (|5.2(l implies (|5.3|l follows from H5.1(l by taking m = 1 and integra- 
ting over the simplex S^. Arguing as in the proof of Corollary 13 one can see that 
it is actually enough to prove (|5.2I) for m = k and k G {1, . . . ,rt — 1}. Note also 
that in the case of equal charges Conjecture ^ would provide a whole new series 
of inequalities for averages of generalized inertia moments involving the critical 
points and the zeros of an arbitrary complex polynomial. These would substantially 
improve both Corollary [3 and the inequalities conjectured in |3D| Remark 2.3.9]. 

Remark 22. In view of the proofs of j2()l Theorem 4.7] and Theorem El above, a 
natural way of attacking Conjecture ^ would be to interpret the arguments of $ 
in (|5.2|) as eigenvalues of some appropriately defined multilinear power "a'°A" of 
a normal matrix A and one of its principal submatrices, where "A" is a suitable 
multilinear functional on ®^Ti. However, defining such a functional appears to be 
a non-trivial question. Note for instance that an obvious polarization of the usual 
antisymmetric tensor power, i.e., setting 



AAB = - 
2 



(A + B)(2)-A(2)-B 



(2) 



for square matrices A, B of the same order, cannot lead to the desired conclusion 
since to begin with one does not get the correct dimensions (these should agree 
with the number of fc-tuples in the two sides of (|5.2|l ). 
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6. Zeros of Borel series with positive 



-COEFFICIENTS 



Theorem 13 applies essentially to any bounded discrete set of positively charged 
particles with finite total charge for which equilibrium points are known to exist. 
We conjecture that if the particles accumulate at a point on the boundary of a disk 
containing all the particles then there exist infinitely many equilibrium points. 

Conjecture 2. Let {ai}i^fi be a sequence of positive numbers whose sum is finite 
and {zijigN a sequence of distinct complex numbers satisfying 



has an infinite set of zeros in the disk B(p) = {z G C : |z| < p}. 

Note that the assumptions of Conjecture|21imply that condition H3.14() holds and 
thus / is meromorphic in B(p). Functions representable as a series of the form 
are known as Borel series and so Conjecture |21 may be rephrased as follows: 
every Borel series with positive l^- coefficients whose poles lie in an open disk and 
accumulate at a point on the boundary of this disk has infinitely many zeros. 

Conjecture El may also be translated in operator theoretic terms by using the 
methods employed in the proofs of Theorems Eld Note first that up to a similarity 
transformation of the complex plane we may assume that ^ = p = 1 . Then one can 
show that the above conjecture is equivalent to the following statement. 

Conjecture 3. Let {zn}nefi a sequence of distinct complex numbers in the 
open unit disk converging to 1 and denote by A the (bounded) normal operator on 
a separable infinite- dimensional complex Hilbert space Ti. satisfying Aen = ZnSn, 
n e N, where {e„}„gN is an orthonormal basis of Ti. If P is a bounded operator 
on TL such that I ~ P is a rank one orthoprojection then the compression PAP\p-h 
has infinite point spectrum. 

Remark 23. The hypotheses in ConjectureOimply that the normal operator I— A is 
compact, which may facilitate its study. However, it seems likely that Conjecture^ 
should actually hold even under the less restrictive conditions \zi\ < p for alH € N 
and linii^^oo = p. In other words, it is reasonable to believe that the convergence 
condition in Conjecture (21 may be relaxed as follows: |z„| < 1 for all n e N and 
lim„_oo \Zn\^l. 

7. Operator versions of the Clunie-Eremenko-Rossi conjecture 

Various results concerning the existence of equilibrium points for Newtonian and 
logarithmic potentials have been obtained in e.g. [HI, 52 ^^'^ E2- I* should be em- 
phasized though that these deal almost exclusively with unbounded discrete charge 
configurations. Note that Conjecture (21 may actually be viewed as a natural analog 
for bounded discrete charge configurations of the following well-known conjecture of 
Clunie-Eremenko-Rossi, see, e.g., [HI Conjecture 2.7] and [5^, Conjecture 1.1]. 

Conjecture 4 (Clunie-Eremenko-Rossi conjecture). Let {aijigN be a sequence of 
positive numbers and {zi}i^fi a sequence of distinct complex numbers such that 



\zi\ < p for i G N and lim Zi — ^ 



^ with 1^1 = p 



for some fixed positive number p. Then the function 



oo 




(6.1) 
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Then the meromorphic function 

oo 

has infinitely many zeros. 

A catchy albeit somewhat loose rephrasing of this conjecture is as follows: every 
flat universe has infinitely many resting points. 

Let us use the framework of i|4.3l to give an operator theoretic interpretation and 
generalization of Conjecture 01 Note first that up to a translation of the variable z 
and a rescaling of the coefficients a„ we may assume wlog that all z„ are non-zero 
and satisfy the normalization condition 



i=i ' 

As in we use the Hilbert space H — P(N) with standard inner product (•, •) 
and standard complete orthonormal system {e„}„gN- Set z„ = r„e'^", n € N, and 
let A and B be the bounded normal operators on Ti. satisfying 

Ae„ = e'''"e„, Bn^r-^e„, n e N. (7.2) 

Clearly, A is unitary and B is compact. Condition 1)7. 1|1 allows us to define a unit 
vector V G 7i by setting 

oo I 

n— 1 ' ' ' 

and one can easily verify that 

{{A-zB)-\,w) ^~f{z), zeC\{z„}„6N. (7.3) 

Set K, = v-*- and let P be the orthoprojection on /C. We use the notation A' and 
B' for the compressions to /C of A and B, respectively, that is, 

A' - PAP\,c, B' - P^P\k. 

By definition, the set of generalized eigenvalues of A' with respect to B' is 

S,,(A') = {zeC|7V(A'-zS')^{0}}, 

where N{T) denotes as usual the null space of an operator T. Let Z{f) be the zero 
set of /. 

Lemma 9. In the above notation one has Z{f) = E^, (A'). 

Proof. We make an appropriate adaption of the argument used in the proof of 
Lemma 13 Using the fact that the 2„ are distinct one can check that if 7^ u S 
then u e N{A— (B) if and only if C = and u = e„ for some n € N. However 



which shows that N{A — (^B) D K, = {0} for C, E C Let now u e /C be such that 
{A' - CB')u = for some ( eC. From the identity 

{A - C-S)u ^{I~P){A~- CB)u + {A' ~ CB')u ^ {{A - CB)u, v)v 

and {{A-CB)u,v) 7^ it follows that {A-CBy^v = {{A- CB)u,v)-^u. By (|73|l 
one gets 

/(C)=-((A-Ci?)u,v)-i(u,v)=0, 
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SO that ( e Conversely, let w E Z{f). Then w ^ {znjnGN and so A — wB is 

invertible by (|7.2|l . The vector u := (A — wB)~^v is clearly non-zero and one has 
(u, v) = f{w) = hence u G /C. Thus 

V = (A - u)S)u = {{A - ■u;B)u, v)v + {A' - wB')u, 

which shows that {A' - wB')u e n {v} = {0} and therefore w <E S^, {A'). □ 

Recall that a vector with all non-zero coordinates in the standard basis {e„}„gi!} 
of TC is called totally non-zero. The above discussion shows that Clunie-Eremenko- 
Rossi conjecture (Conjecture is in fact equivalent to the following statement. 

Conjecture 5. Let {an}neN o,nd be a sequence of unimodular complex 

numbers and a sequence of positive numbers converging to 0, respectively. Assume 
that a„/3~^ are distinct and define a unitary operator A and a compact selfadjoint 
operator B on Ti. by setting 

Aen = anGn, Bsn = /3„e„, n £N. 

Let v be a totally non-zero unit vector in Ti., denote by P the orthoprojection on 
/C = and let A' = PAP\ic and B' = PBP\k.. Then \Y.^,{A')\ = oo, i.e., there 
exist infinitely many z G C such that N{A' — zB') ^ {0}. 

If the vector v in Conjecture El is represented in the standard orthonormal basis 
of Ti as V = ^n^n then Keldysh's value distribution theorem for meromorphic 

functions implies that ConjectureElis true provided that X^^^i /^n ^knP < ^ {see, 
e.g., and [18, Ch. 5, Theorem 6.2]. 

Conjecture |S1 suggests the following more general question: let A be an injective 
compact normal Hilbert space operator and let B and C be finite-rank operators. 
Is it true that |S]j^^c(A + = oo, that is, are there infinitely many z G C such 
that N {A-\- B — z{I + C)) ^ {0}? It is not difficult to see that an affirmative answer 
to this question would actually imply the validity of the operator theoretic version 
of the Clunie-Eremenko- Rossi conjecture (Conjecture [S}. Indeed, in the notation 
of Conjecture [S] we may write 

A' - zB' = -zA{A-^B - A-^[PB{L - P) - [I - P)B] 

- [/ - A-\PA{I -P) + (L- P)A)] } . 

It remains to note that under the assumptions of Conjecture El the operator A^^B 
is injective, compact and normal while both A^^[PB{I — P) — {L — P)B] and 
A^^[PA{L — P) + {I — P)A] have rank at most two. However, the answer to the 
question raised above is negative, as one can see from the following example. 

Example 1. Let J be the simple dissipative Volterra operator on L^{0, 1) defined 
by 

iJf)ix)^2i f f{t)dt 



and set A = J+J*. B = J ~ J* and C — Q. Then A is injective compact selfadjoint, 
B is a one-dimensional operator while the point spectrum of ^ -I- S = 2J is empty 
(cf. ^1 p. 187]), so that in this case A + B — z{I -\- C) is injective for all z G C. 

Let us finally mention that the existence of equilibrium points for some discrete 
potentials generated by complex charges was established in e.g. [22 . We propose 
the following analog of Conjecture 0] for certain potentials of this type. 

Conjecture 6. Let {a„}„gN and {z„}„gN be sequences of complex numbers such 
that the z„ are distinct and satisfy 



oo 

0-nZn > for 71 G N, Z„ — > OO OS 71 — > OO , 

n=l 



< OO. 
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Then the meromorphic function 



oo 



'TL 



z ^ z. 



'n 



has infinitely many zeros. 

As in the case of Conjectures|21and21 one can reduce the conjecture stated above 
to a purely operator theoretic question. Indeed, in the notations of Conjecture 
let A denote the injective compact normal operator on 7i satisfying ^e„ = z,7^e„, 
n € N. Define a vector 



Arguing as in the proof of Lemma IHl one can show that e 2.{f) for any C, € 
Ep(A')\{0}, where Z{f) is the zero set of / and Yjp{A') denotes the point spectrum 
of A' . For Conjecture ini it would therefore be enough to show that |I]p(A')| = oo. 



The author is grateful to Mihai Putinar, Nikolai Nikolski, Alexandre Eremenko, 
Donald Sarason and Boris Shapiro for their interest in this work and valuable 
remarks as well as to Semen Malamud for making available his inspiring paper |26| . 
Thanks are also due to the anonymous referee for numerous constructive suggestions 
and for pointing out several useful references. 
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